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In this paper, I show that backward (in time) running states can be explicitly accounted for by expanding the 
interpolator basis in the variational method in lattice QCD. The backward running states can then be removed by 
choosing an appropriate linear combination of interpolators, which improves the signal significantly. The proof 
of principle, which also makes use of the Time-Shift Trick (Generalized Pencil-of-Functions method), will be 
delivered at an example on a 64'* lattice close to the physical pion mass. 
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Introduction. — In lattice QCD, the variational method 
DM has become a widely used method to separate physical 
states. Its main application is spectroscopy (recent examples 
are, e.g., sni), but it has also been successfully used to ex¬ 
tract specific physical states in hadron stmcture calculations 
(e.g., DSDll). In this paper, I show that backward (in time) 
ranning states can be explicitly accounted for by expanding 
the interpolator basis in the variational method. This is impor¬ 
tant for smdies where the backward propagating states are not 
negligible compared to the forward propagating ones in the 
time range of interest. The backward propagating states can 
then be eliminated with a properly chosen linear combination 
of interpolators. For the proof of principle, the additional in¬ 
terpolators are obtained “for free” using the Time-Shift Trick 
(Generalized Pencil-of-Function method), which I will also 
briefly review in this paper. 

The Variational Method. — The standard variational 
method uses several interpolators Xi,i = 1,..., Nx to sep¬ 
arate the physical states Si,i = 1,..., iV^, with energies Ei. 
Typically, one assumes that the number of significantly con¬ 
tributing physical states is equal to the number of interpola¬ 
tors Nx = Ng = N. The interpolators couple to the physical 
states through 

(s,|2j|0) = tty. 

On a lattice with (anti)periodic boundary conditions, the cor¬ 
relation matrix is 

= ^(0|Ji|sfc)(sfc|2j|0) exp{-Ekt) 

k 

+ X(0|^]|sfe)(sfc|X,|0) expi-EkiT - t)) 

k 

= X Aafci exp{-Ekt) 

k 

+ ^4kbkjeM-Ek{T-t)), ( 1 ) 

k 

where T is the temporal extent of the lattice and I have intro¬ 
duced bij = (0|zj|si). The first term in Eq. ([^ contains the 


forward moving and the second term the backward moving 
states. Neglecting the backward moving states (which is often 
a good approximation), Eq. ([T]) becomes 

Cii(i) = X exp{-Ekt). 

k 

One proceeds by solving the generalized eigenvalue problem 

j 3 

where I denote the fc* eigenvalue and -vector with and 

ik) 

, respectively. This is equivalent to finding the eigenvec¬ 
tors and -values of the matrix 

Gij{to,ti) = ''^C-i}(to)Ckj{ti) ( 2 ) 

k 

= X exp(-E;fc(fi - to)), 

k 

One sees immediately that the eigenvalues are exp(—— 
fo)) and the eigenvector, with components j, is = 
a“/. We can then construct optimal interpolators = 

Xj f which couple - at least in principle - only to the 

physical state sp. 

(sfcixipho) = ^(sfc|x]|0)z;j*^ = 

3 3 

While this method works extremely well in many cases (see, 
e.g., lEHISl), it fails for lattices with (anti)periodic boundary 
conditions when the backward running states are so light that 
they do not decay sufficiently fast and contribute signihcantly 
on the left side. In Eigure[2 I show this problem for pions on 
a 64* ensemble with Wilson gauge action, Nf — 2 flavors of 
dynamical Wilson (Clover) fermions, (3 = 5.29 (correspond¬ 
ing to a « 0.07 fm) and k = 0.13640, where the pion has a 
mass of approximately 150 MeV. Eor ensembles like this one, 
but also at much larger pion masses, the variational method 
has to be modified to yield reasonable results. 
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FIG. 1: The mass eigenvalues, found with the traditional variational 
method with two interpolators, applied to the pion. For t > 9, this 
method extracts the forward (green circles) and backward (red stars) 
running pion. For f < 7, it finds the pion and an excited state (blue 
diamonds), but the excited state does not form a plateau. For t = 8, 
the identification of the found states is unclear (black crosses). 


The cosh Method. — The straightforward way to account 
for the backward running states is to replace the exponential 
decays by coshs (and sinhs, see below). As an example, I 
use two interpolators for the pion (tt) and excited pion (tt*) 
states, Ii = If = d'y^u and I 2 = = djblAU. The adjoint 

interpolators, which couple to the antiparticle, are; 

n'^ = {d'jsuy = U'ysd 

= {dl5lAu)^ = U'jA'lbd = -M7574d- 

In simulations in which the up and down quarks are identical 
(i.e., most of the current simulations), one finds 

(7r+|tf|0) = (7r-|n|0) 

(tt+IAJIO) = -(7r-|Ao|0), 

and equivalently for the excited pion. In other words, b\^ = 
an and = — 0 * 2 - The relative minus sign between these ex¬ 
pressions is responsible for the fact that (Ao|n^) and (IflAj) 
have a sinh-like time dependence, while (n|n^^) and (Ao|AJ) 
are cosh-like. Assuming that the Uij are real, Eq. Q becomes 

(n(f) (0)) = 2a?i exp{-E^T/2) cosh{E^t) 

+ 2021 exp(—£'*T/2) cosh(£'*t) 
(Ao(f)|Aj(0)) = 2012 exp(-T;^r/2) cosh(T;^t) 

+ 2022 exp(—/2) cosh(£^*f) 
(Ao(f)ln'f(O)) = - 2oiiOi2exp(-T;^r/2)sinh(£'^f) 

— 2021 O 22 exp{—E^T/2) sinh{E^t) 

= (n(f)|4(o)), (3) 


the solution is numerically not very stable and needs some 
“supervision”. Once the o^ have been found, one can again 
construct optimal interpolators for the forward ranning states 
by finds for an arbitrary operator O: 

(O(f)|I«^(0)) 

= ^(O(f)|x](0))o-/ 

3 

= ^(0|C>|sfe)(sfe|X]|0)oj/ e-xp^-Ekt) 
k,j 

+ ^(0|x]|sfc)(sfc|O|0)o-/ exp(-T;fc(r - t)) 

k,3 

= {0\O\s,) exp{-Eit) 

+ '^{sk\O\0)bkjaJ^ exp{-Ek{T - t)), 
k,j 


and the example above showed that ba~^ is not necessarily 
the unit matrix. Therefore, while the forward running contri¬ 
bution is optimized for the wanted state, the backward running 
part contains several or all states, which prevents simple ex¬ 
ponential or cosh fits. 

The Variational Method with an Expanded Interpolator 
Basis .— An improvement over the cosh method can be 
achieved by using the standard variational method with an ex¬ 
panded interpolator basis, where the additional interpolators 
account for the backward running states. In other words, I 
treat the backward running states as forward running particles 
with negative mass. Usually, the heavier backward running 
states will be negligible on the left side since they decay too 
fast. The number of backward running states that contribute 
significantly is Ni, < Ng. One then needs Nx = Ng + Nf, 
interpolators to account for the forward and the significant 
backward running states. I define 6^^ = b^j exp{—EkT/2) 
and sort the states so that the first Ni, states are the ones that 
are non-negligible to rewrite Eq. Q as 


Ns 

exp{-Ekt) 

k^l 

iVb ^ ^ 

+ exp{Ekt). 

fc=i 

With 

, ^ f o,kj, 1 < k< Ns 

"^■“1 Ns<k<Nx 

and 

p ^ l<k<Ns 

^'^-{-Ek-Nn Ns<k<Nx, 


where I have used E^ = E^, and t = t - r/2 for brevity. this becomes the standard variational problem. I define the 
Eqs. ([^ have to be solved numerically and it turns out that optimal interpolator as Then, for an arbi- 
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trary operator O, one has (dropping the negligible states) and therefore 


Ns 

« g ^(0|O|sfe)(sfc|z]|0)d-/ eM-Ekt) 


k=l j 
Nh 

+ ^^(0|z]|sfc)(sfc|O|0)d-/exp(-£;fe(T-f)) 
k^l j 
Ns 

= y^(O|0|sfc)afcjrf~/ exp(-£'fct) 
k=l j 
Nh 

+ y^(sfc|O|0) exp{-EkT/2)bkjdJi^ exp{Ekt) 

k=l j 

(0|C>|si) exp{-Eit), I <i< Ns 

(si_ArJC>|0) exp(Ei_N,i^, Ns < i< Nx- 


This shows that the variational method with an expanded in¬ 
terpolator basis works indeed and one is left with an optimal 
interpolator for the ground or excited states, with no contam¬ 
ination from the backward running states, which greatly im¬ 
proves the further analysis. 

In my example, I need Nx = 3 interpolators to account 
for the pion, the excited pion and the backward running pion. 
These interpolators can be obtained using different currents, 
various levels of smearing or any other method that yields lin¬ 
early independent interpolators. Here, only data for the two 
interpolators H and Aq with one level of smearing each was 
available. So in order to avoid using additional supercomputer 
time, I had to use the Time-Shift Trick to obtain the third in¬ 
terpolator. 

The Time-Shift Trick .— The Generalized Pencil-of- 
Function method ll20l or, more catchy, the Time-Shift Trick 
allows one to obtain additional interpolators “for free”. As 
the name suggests, new interpolators are constructed by time- 
shifting other interpolators. I will label the time-shift operator 
T{5t) = T^St), with the amount of time shift 6t. Its action 
on an operator 0{t) is defined by T{St)0{t) = 0(t —St) and 
due to time invariance one sees that 


{02{t2)T{dt)0i{h)) ={02{t2)Oi{ti-St)} 
= {02{t2+St)0i{h)) 


0{t)T{St) = 0{t + dt). (4) 

Using this, one can define a new interpolator I' = T{St)Xi 
and sees immediately that 


{s,\I'\0) ={s,\Ti6t)IfO) =exp{-E,6t){s,\IfO) 

= exp{—EiSt)aij. 

Therefore, I' is, in general, a new linearly independent inter¬ 
polator. 
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FIG. 2: Comparison of the cosh method (filled symbols, slightly 
shifted to the right) with the variational method with an expanded 
interpolator basis (open symbols). The backward running pion, the 
forward running pion and the excited pion state are denoted with red 
stars, green circles and blue diamonds, respectively. We note that 
both methods find the pion and the first excited state, even with simi¬ 
lar error bars. The variational method with an expanded interpolator 
basis is numerically more robust than the cosh method and therefore 
able to go two time steps further to the right. 
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Proof of principle .— In this final paragraph, I show that 
the variational method with an expanded interpolator basis 
and the Time-Shift Trick actually work with “real data”. I in¬ 
vestigate the pion with the two interpolators H and Aq and the 
time-shifted interpolator 11 T{6t) so that I can account for the 
pion, the excited pion and the backward running pion. Using 
Eq. (0, one obtains the correlation matrix 


/ (n(f)int(o)) (n(f)| 4 ( 0 )) {u{t + st)\u^m \ 

C(f)= (Ao(f)int(O)) (Ao(f)| 4 ( 0 )) (Ao(f + <5<)|nt(0)) . 

V (n(t + 5f)|nt(o)) (n(f + ^f)|4(o)) {u{t + 2 St)\u^o)) J 


I use this matrix in the standard variational method and find suits, I vary to, ti and St. It turns out that for this example, 
the eigenvalues of G{to, ti) in Eq. Q. To obtain optimal re- ti—tg = 2 and St = 14 are good choices. Eor the sake of com- 
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FIG. 3: The correlators of the optimal interpolators with the 11 in¬ 
terpolator, based on the cosh method (filled symbols) and the AIM 
(open symbols). The lines are best linear fit lines to the log of the 
correlators obtained with the variational method with the expanded 
interpolator basis in the region where the excited states have decayed 
sufficiently. The backward running pion (red stars) is by construction 
only found with the variational method with an expanded interpola¬ 
tor basis. The forward running pion (green circles) is influenced by 
the backward running state for the cosh method while it is a per¬ 
fect exponential up to large times for the variational method with 
an expanded interpolator basis. Also, the excited pion state (blue 
diamonds) is resolved better with the variational method with an ex¬ 
panded interpolator basis than with the cosh method. 


pleteness, I compare these results to the results obtained with 
the cosh method, where I also kept ti — Iq = 2. The results 
are shown in Figure Both methods find perfect plateaux for 
the pion (the variational method with an expanded interpola¬ 
tor basis also for the backward running pion) and reasonably 
nice plateaux for the excited pion. The error bars are similar 
and the eigenvalues agree within the error bars in the region 
of the plateau. So while the results are comparable, I prefer 
to use the variational method with an expanded interpolator 
basis over the cosh method because of its numerical stability 
and robustness (which also allows the former method to reach 
higher to than the cosh method). Furthermore, the variational 
method with an expanded interpolator basis is faster than the 
cosh method, although this is usually not relevant with today’s 
computers. 

The big advantage of the variational method with an ex¬ 
panded interpolator basis over the cosh method, however, are 
the optimal interpolators. For each method, I construct the 
optimal interpolators as described above and show the correla¬ 
tors (n(f)|lo^pj^(0)) in Figure]^ The variational method with 
an expanded interpolator basis yields an optimal pion, i.e., one 
that has an exponential (linear in the logarithmic plot) behav¬ 
ior throughout the plot range. The backward running pion is 
also very good and the excited pion is reasonable. The cosh 
method, on the other hand, yields a poor excited state and a 
pion which is usable up to f « 10. For larger t, the backward 
running state starts to strongly influence the signal. 


In conclusion, the variational method with an expanded in¬ 
terpolator basis has proven to be a robust method that is at 
least as good as the cosh method in finding the energy eigen¬ 
values but has its real strengths in the construction of optimal 
interpolators. The Time-Shift Trick is also very useful and 
can be applied wherever additional interpolators are needed, 
in many cases even at no extra computational cost. 
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